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Abstract 

We study the behavior of the renormalized sextic coupling at the intermediate and strong coupling 
regime for the theory defined in d = 2-dimension. We found a good agreement with the results 
obtained by the field-theoretical renormalization-group in the Ising limit. In this work we use the 
lattice regularization method. 


1 Introduction 


The scalar field {>^'^)d =2 theory is known to have an interactive continuum limit. This theory is 
an excellent theoretical laboratory in which nonperturbative methods can be tested. For theo¬ 
ries in less than four space-time dimensions can offer interesting and complex behavior as well as 
tractability, and for example the case of three space-time dimensions, they can even be directly 
physical, describing various planar condensed matter systems. The 99 ^ theory has been study in 
the regimen of the intermediate and strong coupling analytically using the path integral represen¬ 
tation of the generating functional of the Green’s function. Two of these methods are the high 
temperature series approach of Baker and Kincaid [1] and the strong coupling expansion of Bender 
et al [2, 3]. Each of these methods is based in the replacing the continuum by a Euclidean lattice 
and expanding the kinetic energy. For example in the strong coupling expansion the method con¬ 
sists in expanding the Green’s functions in the inverse powers of the bare coupling. Of course the 
result of this expansion will be better as the bare coupling is strong enough. This is because the 
expansion requires that the coupling constant has to be large to perform the analytic expansion of 
the path integral. Another nonperturbative analytic method which was proved to be a powerful 
tool is the the field theoretical renormalization group [4]. 

The first study of the 99 ^ theory in the regimen of strong coupling using lattice Monte Garlo 
simulation was done by Gooper, Freedman and Preston [5]. They showed that the theory has a 
non zero scattering amplitude in one and two dimensions and the asymptotic value for the renor¬ 
malized coupling constant was in good agreement with the high temperature series and strong 
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coupling expansion method. Freedman et al [6] continued the study in three and four dimensions 
and they found Monte Carlo evidence that the lattice theory has a noninteracting contin¬ 

uum limit consistent with the result of renormalization group calculation [4] and strong coupling 
expansion [2, 3]. In contrast the {ip^)d =3 theory appears to have an interacting continunm limit. 
In this paper we are interested in the higher order renormalized coupling constant in particular 
the sextic coupling constant g^. We expect that in the limit of strong conpling this quantity will 
approach the asymptotic limit. This because the theory in the strong conpling limit can be consid¬ 
ered to be in the the critical region where fluctuations are so strong that they complete screen ont 
the initial interaction so that the behavior of the system becomes universal. In two dimensions, 
Sokolov and Orlov used renormalization group expansion and Pade-Borel-Leroy resnmation tech¬ 
nique to get qq [8]. In the literatnre there are few works available of the sextic coupling constant 
using the approach of lattice Monte Carlo. The hrst work was performed by Wheter [9] and later 
by Tsypin [10]. So it seems appropriate the nonperturbative study using the lattice Monte Carlo 
techniqne [11] to get not only the asymptotic valne of the higher order renormalized conpling 
constant but also the behavior on the intermediate and strong regimen. Here as usual we consider 
h = c = 1. 
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2 The effective potential 


We consider here the 0^ theory in d-dimensions Euclidean space in the presence of a source J, 
which the bare action is given by: 


J] = J d^x 


+ ^00 ~ "^000 


We introduce the vacuum persistent function Z[J]: 


Z[J] = J 


( 1 ) 


and from this we obtain the connected n-point Green’s functions Gn(xi.. .Xr. 


. . . Xyi) 
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5J{xi) 5J{Xn) 


\nZ{J) 


( 2 ) 


j=o 


The renormalized proper iV-point Green’s function is obtained as follows: the wave- 

function renormalization is obtained from the two-point Green’s function from: 


= 


and the renormalized mass is dehned by 




dp^ 


( 3 ) 


p2=o 


m; = ZG 2 (p ) 


p2=o 


( 4 ) 


The effective action r[0c] is dehned by a functional Legendre transform of W[J] = lnZ[J]: 


r[0,] = I d^x 0c(a;) J(a;) - W[J], 


( 5 ) 
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It is well known that the effective action is the generator of proper Green’s functions and in 


particular we can write the effective action in the form 

1 


r[0c] f . . .,Xn)4>c{Xl) ■ ..4>c{Xn)- 

„ n J 


( 6 ) 


Here , Xn) is the proper n-point Green’s functions in presence of the source J. We 

dehne the effective potential, G(0), by 


rui = / d^x 




(7) 


an the renormalized effective potential 




r<2"i(o) 


i2n 


2n! - 

where = (j) and the renormalized r’^”^(0) proper n-point Green’s functions are given by 


( 8 ) 


rW(0) = 


(9) 


from here we see that Ur is the generating function of one particle irreducible renormalized Green’s 
function at zero external momentum on all legs. 

The particular interest to us are the renormalized coupling constant r(^)(0) and the renormalized 
sextic coupling constant f(®)(0) which can be expressed as follows: 


rW(o) 




p=0 


and 


rf (0) = -Z^G^^f(Ge - lOGlG^') 


( 10 ) 


( 11 ) 
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The quantity that will be extracted from lattice Monte Carlo simulation is the dimensionless 
renormalized constant g 2 n defined by, 


g2n — 


r(2-)(o) 


ml^-nd+d(^2n) 


For d = 2 this expression becomes 

mj{2n)\ 

and the effective potential divided by ml can be written, 


92n — 


Ur{(j)r)/ml = -(j)l + g4(j)t + + 


( 12 ) 


(13) 


(14) 


3 Simulations results 

For the Monte Carlo simulation of the lattice field theory described by Eq. (2) we used the standard 
Metropolis algorithm. In order to avoid the trapping into meta stable states due to the underlying 
Ising dynamics we use the cluster algorithm using the embedded dynamics for 0^ theory, according 
to the action, [7]. 

*S*Ising -EE \(j){x + s(x + fl)s{x) , (15) 

a: fi 

where s{x) = sign(0(x)). Statistical errors are evaluated taking into account the autocorrelation 
time in the statistical sample generated by the Monte Carlo simulation. In the simulation we used 
the lattices 16^, 20^, 36^ and 64^. To minimize finite size effects we choose 1 << l/(amr.) << L. 
For the purpose of numerical simulation it is appropriated to work with equivalence Lagrange on 
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the lattice 


•SM = 2(t ^ 4>{x)4>{x + e^) + ^ Ifi(x)^ + Yi MMxf - 1)' 


(16) 


X,fl 


The parameters and are related to k and A as follows, 


2 1 - 2A 6 A 

m = - 2d g = 

K K,^ 


(17) 


We fonnd the bar parameters {m^,g) for the different lattices by hxing the valne of A and varying 
the parameter k until we hnd an acceptable value of the renormalized mass within few percent. 
All the values of bare mass were negative. Of course we can not have a phase transition since we 
imposed the renormalized mass to the condition 1 << l/{amr) « L. The renormalized mass 
nir was calculated using the expression 


-2 -2 
a rrij. = 


(L\" 

r<0(o)2>-< i0(p)|2>i 

UttJ 

< |0(P)P > 


(18) 


where here (p is the Fourier transform of (p and p = (x^,0) is the smallest available non-zero 
momentum. 

For purpose of calculating the values of 5^4 and g^ from simulation we use lattice version of 
eqs.(lO) and (11), 


^ 4 ( 4 !) = 


FW(0) 


;< 0(0)^ > -3 < 0(0)^ 

< 0 ( 0)2 >2 


= —m: 


mt 


(19) 


^6(6!) = 


F(6)(0) 


mt 




mt 


^^2^2 <<P{^f > -15 < 0(0)" >< 0(0)^ > +30 < 0(0)^ >3 


< 0 ( 0)2 >3 
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L 


9& 


9a 

16 0.60 ±0.013 1.18 ±0.06 

20 0.60 ±0.011 1.16 ±0.04 

36 0.59 ±0.012 1.12 ±0.05 

64 0.585 ±0.012 1.112 ±0.04 

Table 1: 

where < 0(0) > is the Fourier transform evaluated at zero momentum . We notice that eqs.(19) 
and ( 20 ) are efficient in the intermediate and strong conpling regime where the statistical errors are 
reasonable contrarily to the case of weak coupling. In table 1. we present the results of simulations 
for various lattices size for a limit of strong coupling. The results in table 1. are consistent with 
the values obtained by the held-theoretical renormalization-group in the Ising limit 5^4 = 0.6125 
and qq = 1.10 [ 8 ]. In Fig. 1. we present the resnlts of simnlation for and for a large range of 
the bare coupling constant (gf ~ 15 — 3232). We observe both g^ and g^ approach to asymptotic 
constant value as we expected. From the hgure also we can see that the variation of both qnantities 
are small considering the large variation on g. 

4 Conclusion 

In this paper we analyzed the behavior of the renormalized sextic conpling for the theory ( 0 ^)rf =2 
on the lattice at strong conpling constant. We fonnd that this qnantity tend to constant value 
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Figure 1: The behavior of and qq with the bare coupling constant 


that is in good agreement with value obtained by using the field-theoretical renormalization-group 
in the Ising limit at four loops. We also notice that the qualitative behavior of is similar to 5 ^ 4 . 
We expect that this behavior is also similar for the case d = 3, which will be presented on future 



report. 
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